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ABSTRACT 


A procedure is presented for solving the problem of allocating air- 
strikes for the interdiction of a transportation network when the flow of 
supplies is limited by the capacity of the system. The damage function is 
assumed to be exponential. 

The inputs required are the upper and lower bounds on the capacity of 
each arc, the vulnerability of the arc to attack and the number of air- 
craft sorties available for the mission. The procedure determines the 
segments in the network to attack as well as the level of attack to re- 


duce the network flow capacity to a minimum. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


an arc between nodes i and j 

maximum length (capacity in the primal) of arc (i,j) 
a lower bound on the length (capacity) of arc (i,j) 
the vulnerable length (capacity) 


the number of sorties allocated to arc (i,j), k 0 


‘ne 
ij - 
the resource constraint, or total units of effort 
available 


a measure of the vulnerability of arc (i,j), or a 
measure of the efficiency of the weapon system used to 
attack arc (i,j) 


the surviving portion of the vulnerable length 


the resulting length of arc (i,j) when kK. sorties are 
allocated. J 


the total number of arcs in a route from source to sink 
the total number of nodes in the network; the sink 
actual flow in arc (i,j) in the primal 

the Lagrange multiplier 


>: (1 Abyss es) eee 11) 
W. Db..-U "J JM 
fe 


the value of u such that h(u) = k 


h 


the qt largest value of W..b:.. on a route 


Ue) 
the minimum feasible length of the route being considered 


proportional damage parameter 


the number of distinct routes from source to sink without 
loops 





IT. INTRODUCTION 


A. SUPPLY NETWORK 

In order to conduct sustained military operations any distance from 
its permanent base a military force must have some means of resupply. 
Quite often a major portion of this resupply comes over the land trans- 
portation system. This is particularly the case when the opposing forces 
may have air and/or sea supremacy or the means of resupply by sea or air 
are not available. 

The ability of the military force to operate at various levels of 
combat performance may be seriously affected by the capability of the 
transportation system to carry the necessary supplies. Limitations on 
the cargo-carrying capacity of the transportation network may occur as 
the result of the characteristics of the road network and vehicles being 
used.. This capacity, expressable in terms of tons per day, measures the 
ability of communications lines to transport cargo. 

The type and amount of supplies required to sustain a military force 
vary greatly depending on the type of warfare being conducted as well as 
the level of combat. Small guerrilla forces do not require as much as 
larger conventional forces and may be able to provide much of their own 
Support from the area they are working in. In general however there are 
military supplies which must be brought in. Logistics planning must 
consider the supply capacities of roads and transportations requirements 
to support a military operation. 

The Army Field Manual for Staff Officers on Organization, Technical 
and Logistical Data [1] lists and describes planning factors for the logis- 


tical support of various size forces in different types of environment. 


The requirements to support operations in various situations are ex- 
pressed in terms of pounds-per-man-per-day and include all classes of 
Supply. 

Holliday [2] has provided the research allowing the concept of net- 
work capacity to be applied to highways. In work done for the Advanced 
Research Projects Agency, he provides the methodology for determining 
road capacities. He analyzes traffic cargo flow based on operational 
factors and relates the road supply capacity to the number of vehicles 


required to achieve that capacity. 


B. INTERDICTION 

When the resupply capability of one force is limited by the capacity 
of the transportation network over which it must operate, it is to the 
advantage of the opposing force to consider the interdiction of that net- 
work. It may be worthwhile diverting resources toward network interdic- 
tion in order to further reduce the capacity of the transportation system. 
By restricting the capacity, the flow of goods is ultimately restricted. 
This is taken as the objective rather than minimizing flow directly. 
Moving vehicles become difficult to interdict when provided with suffi- 
cient defense, jungle canopy, camouflage or cover of darkness. In 
addition, destroying vehicles with cargo may not accomplish the desired 
objective unless either vehicles or supplies are in short supply. 

Airstrikes may be used to further restrict the capacities of segments 
in the network. This may be done by destroying bridges, ferries or rail- 
road track, cratering or blocking roads, et cetera. In addition, aircraft 
may reduce the capacity flow of traffic over a network by reducing the 


period during which vehicles may operate. 


rr — — 


Several authors have considered the supply network interdiction 
problem. Wollmer [3,4] determines the most vital links of a network and 
examines the effect on capacity of removing the most vital link from the 
network. Arcs are subject to breakdown which results in a decrease in 
capacity (1) by a known quantity and then (2) by a quantity which is a 
random variable. He finds the Seaiined reduction of flow possible if 
n breakdowns occur and where these must occur. No consideration is made 
for the cost of the breakdown which, in the case of interdiction, corre- 
sponds to the number of aircraft allocated. In addition his method 
appears cumbersome when there may be more than a few breakdowns per arc. 

Durbin [5] evaluates the ability of a transportation network to 
deliver supplies as the road segments are successively destroyed and 
repaired. The methodology is provided for successively removing the next 
most vital arc until the flow is reduced to zero or the pre-determined 
number of links have been destroyed. The process then steps to the next 
day, restores repaired links and continues. He too uses lower bounds of 
zero and does not consider the cost to destroy the link. 

The effectiveness of airstrikes against a segment of the network may 
be a function of many factors. The type of weapon system involved, its 
accuracy and destructiveness, characteristics of the segment such as 
overhead cover and concealment, defense, vulnerability and recoverability 
are all examples of such factors. 

Mustin [6] has looked at the interdiction problem when the decrease 
in capacity of an arc has a fixed deterministic rate and is known; the 
reduction in capacity is assumed linear between upper and lower bounds 
with the latter being positive. He proposes a linear approximation when 


the returns are not constant. 


1] 


Generally airstrikes can be expected to reduce the capacity of a seqg- 
ment of the network according to some sort of law of decreasing marginal 
returns. This is reasonable when it is considered that although a bridge 
may be destroyed, it may still be possible to ford or ferry supplies 
across a river. Similarly shuttling may be conducted to overcome the 
effect of cratered or blocked roads. 

This paper therefore considers a damage function which exhibits 


decreasing marqinal returns throughout; namely, an exponential function. 


C. ALLOCATION 
The general interdiction problem facing the strike planner on a 
given day 1s which segments of the network to attack and at what level. 
His objective 1s to create the maximum reduction in flow capacity in 
order to have the greatest effect on the opposing side's combat perform- 
ance. The restriction placed upon him is the number of sorties available. 
To accomplish his objective the strike planner must make the most 
efficient utilization of the resources available to him. He must consider 
not only the capacity of each segment of the network but also its vulner- 
ability and importance. An arc with the largest capacity is not necessarily 
the most worthwhile target, nor is an arc with the highest vulnerability. 
It is assumed that the strike planner has available to him the neces- 
Sary information concerning the factors involved. He Knows the upper and 
lower bounds on capacity of each segment and its vulnerability. The number 
of sorties available to him is fixed. His problem is how to allocate those 
sorties in order to have the greatest effect on the enemy's resupply capa- 


bility. 








IT. QBJECTIVE AND SCOPE 


The purpose of this paper is to present a method of solving the 
problem of allocation of effort in the interdiction of a transportation 
network for a damage function having decreasing marginal returns. 

The strike planner is assumed to have available to him information 
regarding the upper and lower bounds of segment capacities. It is also 
assumed that he knows the vulnerabilities of each segment. Finally, the 
total number of sorties available is assumed known and fixed. A necessary 
assumption is that the various segments are independent in so far as the 
reduction of capacity of one arc does not affect the vulnerability of 
another arc. The proportional reduction in capacity of arc (i,j) with 
Ka sorties is assumed to be 1 - exp(-b. sk; 5). 

The solution procedure selects the segments to be attacked which will 
give the smallest resulting value of maximum capacity and the number of 


sorties to allocate to each segment. Also given in the results is the 


minimum value of capacity obtainable with the number of sorties available. 
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A. NETWORK DESCRIPTION 

The transportation system can be represented by a network of arcs 
and nodes. Nodes represent intersections of road segments. Further, 
they may be used to represent any point at which it is convenient to 
distinguish between the road characteristics on either side of the node. 
Arcs represent road segments. Each arc joins two nodes. Associated 
with each arc are quantities representing the upper and lower bounds on 
capacity and the vulnerability of the arc. Both upper and lower bounds 
are in terms of the maximum amount of flow which may pass over the arc 
in a unit of time, say in tons per day. 

Two special nodes are the source and sink. It is assumed that al] 
goods flow from the source to the sink. If there are several sources 
or sinks, this may be taken care of by adding artificial nodes and arcs. 
The capacities of these artificial arcs are infinite and they are not 
vulnerable to attack. 

The notation (i,j) represents the arc between nodes i and j. Nodes 
are numbered from 1 to N with | the source and N the sink. The inter- 
mediate nodes have any of the values between 1 and N. For a particular 
route, n will represent the total] number of arcs in the route. 

The actual flow in an arc (i,j) 1s designated as XG if it is flowing 
from 1 to j and as Xs 4 biol 1s flowing from j to i. Flow is assumed 
to be from the source to the _" but may be in either direction over the 
intermediate arcs and nodes. Due mm the undirected structure of the 
network considered, capacities on arcs represent bounds on net flow over 


the arcs in either direction. 


The principle of flow conservation at the nodes is assumed. This 
means that flow out of a node equals flow into a node; that is, there is 
no storage at nodes. 

The capacity on arc (i,j) at any time is represented by mM; and 
implies 


0 < ie <M, 4° 


As. long as the capacity of an arc is the same in both directions, mi 
=. My, This will be the nature of network arcs considered. 

The lettering Us and Ls represent the upper and lower bounds 
respectively on the capacity of the arc (i,j). If there is no lower 


bound then Li =. 0. The following relationship holds: 


The vulnerability -.of. arc (i,j) is represented by the parameter Ds 3: 
This vulnerability. parameter indicates the efficiency of the particular 
weapon system against. the arc. 


The amount of capacity of arc (i,j) which is vulnerable to attack 


18 Wa = Ue Lig. As the expected proportion of reduction is 
1 - exp.(-b; ky 4), the reduction function which results is then 
W, jl - exp(-b, sk. 5) J 
Was would be the expected reduction in capacity for a very high 


effort; tlnat-is, for.k,. =o. Only by allocating an infinite (or very 
high_amot ant of) effort toward arc (i,j) would the expected resulting 
capacity: be its lower bound. 

For finite kas the remaining capacity of arc (i,j) is given by: 


Mia = U; - Well - exp (-by sky 


ie 


B. MAXIMAL FLOW DETERMINATION 

It is assumed that the opposing force wishes to maximize the total 
flow, Q, through the transportation network available. This is the 
Standard maximal flow problem which can be stated as the following 


linear programming problem: 


maximize Q 


SO = ae 
k lk j JN 
Subject to 2 Xai ; Xj = 0 | ne A , N-1 
ae eee. yN 
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where Mm; = 0 if the particular arc (i,j) does not exist; otherwise ms 
is a known constant. The notation Xap represents the flow from node i 

to any node k. The notation A represents the flow from any node j 
into node i. 

The maximal flow can be determined by use of the maximal-flow 
minimal-cut theorem [7]. A cut set, defined as a set of arcs which 
separates the source from the sink, has a value defined as the sum of 
the capacities of all arcs in the cut set. The maximal-flow minimal-cut 
theorem states that the maximal flow in a network is equal to the minimal 
Value of all cut sets. 

It is convenient to use the notion of a topological dual in finding 
the cut set through a network [8]. The topological dual will be a net- 
work where arcs have lengths instead of capacities. There 1s a one-to- 


one correspondence between the cuts of the original network and the 


routes through the dual. The shortest route through the dual corresponds 
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to the minimal cut set in the original or primal network; its length is 
equal to the value of the minimal-cut set. The topological dual is 
defined only for planar networks, that is, for a network where no two 
arcs intersect except at a node when the network is represented on a 
plane. This will be the only type of network considered in this paper. 

To form the dual, an additional arc is drawn to connect the source 
and sink of the primal. This forms the modified primal, G. The dual 
of the primal is constructed by taking a vertex inside each region of 
G and one vertex outside G and connecting vertices in adjacent regions 
with arcs which cross arcs of the primal. The vertices are points which 
become the nodes of the dual. Only one arc can cross an arc of the primal 
and no arc crosses the artificial arc. Select as the source or sink the 
node outside G and then designate the node in the region created by the 
artificial arc of G as the other end of the network. 

The lengths of the arcs in the dual are the capacities of the arcs 
which they intersect in the primal. Because the primal network is un- 
directed, the arcs of the dual are also undirected. The maximal flow 
through the primal can therefore be determined by finding the shortest 
route through the dual. The shortest route is found by finding the path 
from source to sink such that the summation of arc lengths in the path 


iS minimum among all possible paths from source to sink. 


C. JHE INTERDICTION PROBLEM 

The interdiction problem associated with the dual involves deter- 
mining the feasible distribution of effort over a route which will result 
in reducing its total length by the largest amount. The improvability of 
the arc in the dual is the same as the vulnerability of the primal arc. 


Each unit of effort (each sortie) allocated toward an arc in the dual 


ee 


shortens the length of that arc resulting in the reduction of the asso- 
ciated arc capacity in the primal arc. Some route through the dual will 
have the shortest length after optimal allocation. The problem is to 
find this "shortest" route and determine the optimal allocation to attain 
its final length. Mathematically, the problem is to find that route 
which 


Mimi zes Sere UG - Wg lie exp(-b; sk; 5) J) 


subject to ee eS OR ee 
vnaade 1j- 
T3J 
In the equations above and throughout this paper, the notation 


y 1S taken to mean the summation over all arcs in the particular 
lees) ; 
route. 

The restraint of the problem above can be considered an equality as 
the objective function is continuous and decreasing for finite K; that 
is, the value of the objective function for any number of sorties less 


than K will always be greater than the value of the objective function 


forefitite K. 
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IV. SOLUTION PROCEDURE 


A. PREVIEW 

By utilizing the topological dual, the problem of allocation 
becomes one of determining the shortest feasible route through the dual 
In finding the solution, there are three separate problems which must 
be dealt with: 


1) which route through the dual to select; 
2) which arcs in the route to attack; 
3) how much effort to allocate to each of the attacked arcs. 


The amount of computations is reduced by establishing upper and lower 
bounds on the feasible lengths for each route, thereby eliminating some 
routes from further consideration. 

Although an integer solution is not provided in the calculations of 
the procedure, the value of the objective function obtained in the 
example problem was not changed greatly by using a simple round-off proce- 
dure. This is generally true due to the nature of the exponential 
distribution. In the event that the parameters are such as to be working 
on the steep slope portion of the exponential curve, it may be necessary 


to investigate the integer solutions. 


B. ALLOCATION OF STRIKE EFFORT TO A GIVEN ROUTE 
In approaching the problem, it is useful to consider a problem pro- 
posed by Koopman [9] and explained in more detail by Danskin [10]. 


Consider first the following lemma by Gibbs as given by Danskin: 


Suppose the set x? = ce ae. maximizes 5 F(x.) 


subject to the conditions » x, = X, X, > 0. Suppose further the 
71 
f are differentiable. Then there exists a scalar u, such that 


Fs 


where Fi is the derivative of f with respect ot Xe 
For the allocation of effort toward interdiction the f. of Gibbs’ 
lemma is Wy all - exp(-b, sk. .) J. From the conditions of the lemma, the 
following relations must be satisfied for an optimal allocation on a 
route: 
W, 5d: exp(-b; sk, .) =u if ki, > 03 (1) 


W; sb, <0 Tee oe = 0. (2) 
It will be desirable to show that (1) and (2) are necessary and 


sufficient conditions. Suppose W..b.. < u and ee > 0, then (1) must 


Tadd 
hold. Because exp(-b. «k, «) is always less than 1 for Diy and Ks both 
positive, then WiaP a5 > u, which is a contradiction of the original 


assumptions. Therefore Ks 5 Oe uit Kae > 0 is required for feasibility. 


Thus W..b.. < u implies Kee = 0. This together with (2) implies 


Ie = 
W; 3b, <u if and only if ks =0. (3) 
The feasibility requirement together with (3) indicates that W, 3b. 
> u implies Kaj > 0 as Kee = Q can only happen if W445 ee 
The value of Kegs where ks > 0, may be derived from (1): 
ks = (1/b; 5) In(W, sb; ./u) (4) 
where u must satisfy 
) | . 
ai (1/b; 5) In(W, sb; ,/u) = kK 
"W)C 


The summation sign in the above relation is taken as the summation over 


all arcs in the route with Waa greater than u. 
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To determine the value of u, first define h(u) and * as follows: 


= 7D... aD a 
mu) =z (1/45) Im 3b; 5/4) 
Wj ig 
= 3 A7b..) Wasp...) - 1 rr 
BND, >a ” ig? 7 "J ij? _ '. > ay 
Wek ag 9 
Vo = Wa = the q-th largest value of We gas in the route. 


The desired value u* is that value of u which satisfies h(u) = K. 
It is heloful to look at the nature of the function h(u). For 
u small (near zero), h(u) is large and for u sufficiently smal], h(u) 


is laraer than any finite K. For u greater than or equal to v, (the 


] 
maximum value of W;3b.. on that route), h(u) is equal to zero. For u 
between 0 and Vio h(u) is strictly decreasing and continuous. 


To find u, the following recursive formula is developed for h(va) 


where v <v <v,: 
n—-q- ] 
q-] 
h(v,) = h(vg_y) + Cin(v, 4) - Intv,)] J (1/b.) 
q-2 q-2 
= D. _) - 
where h(v, 4) A a .) In(v.) In(v__4) a (1/b.) 


Specifically, 


h(v_) 


I 
ope. 


* 
” (1/b; ) In(W; 3b; 5/u ) 
igrig 
n = 
Yee ini .o./uU oF (5) 
j=] | 13 
The single subscript on the b's indicates the b, corresponding to the 


Vo which is the q-th largest value of W gb. 5. 


al 


i 
If K > h(v)s then there exists au < Vi Such that 


k= h(v,) + EIn(y fu") 


ims SS 
a | 


(1/b,). 


1 
Then 


* 


u = ov exp{-[K-h(v, ) ]/ (1/b,) } 


e 


2d 
= 


] 


In the event K < h(v)s h(v_) is successively computed for g =n, 


q | 
n-l,...,1 , that 1s, from the smallest to the largest value of ee until 


h(v,) > K ? h(v 4) 


for some q. Then the desired value of u is 


] 


d= 
x (1/b,)) 


EXD {-[K-h(v,_) 2 


The value of Kags the amount of effort to allocate for an arc 


having W..b.. > u, may be determined from (4): | 


i ay 
; In(W. cbs ./U> : 
ka = (1/b; 5) n{W 5 /u Ni, (6): 
* * 
When No ae a EUs SG = 0, as determined in the derivation. © 


It remains to determine which route in a network to select in order 


to get the minimum feasible route through the network. Letting 


Mo o= ¢ {U.. - Ws LT - exp(-b; .k;.) J} 


* 
for a particular route after the Kay have been determined, the problem 


* 
is to find that route for which M is minimum over all possible routes. 


Zee 


For this a stepwise procedure is developed to eliminate some of the 
routes from consideration by means of establishing bounds on the feasible 
route lengths. 

For the special case of - = 0, a simplification results because 
W 


= U The criterion for attack of arc (i,j) is that U; 3b. is 


id 1) 
* 
greater than u. This can be handled by direct substitution and 


presents no problem in the computations. 


C. STEPWISE PROCEDURE 

1) Determine the minimum normal route using the upper bounds on 
lengths, 5 This may be done using any of the shortest route algorithms 
{11}. The route through the network is found such that the sum of the 
upper bounds is the minimum over all routes. This then becomes the first 
least upper bound for the network. 

2) Determine the r-shortest routes through the network using the 
lower bounds on capacities, where r is the number of distinct routes in 
the network from source to sink without loops. 

Because of the nature of the problem it is not necessary to 
consider any route with a lower bound greater than the least upper bound. 
In general, therefore, the number of routes that have to be considered 
will be considerably less than the total number of distinct routes in the 
network. In addition, there is no practical reason for passing over an 
arc or node more than once because, in determining a minimum cut set, if 
vs is the length of a feasible route without loops, any route considering 
loops would have a length greater than Mw and need not be considered. 
Therefore if some computational procedure is used to determine the r-th 
shortest route, caution must be taken that no arc or node is included twice 


in any route. Clarke, Krikorian and Rausen [12] provide a procedure for 


ZS 


determining the N best loopless paths in a network where N is any 
positive integer. 

3) Beginning with the route having the smallest bound, determine 
improved bounds by computing hv ) from equation (5). 

4) As h(va) is a continuous and decreasing function, if K is 
greater than h(v), this means that u is less than 7 and all nt are 


positive. The surviving portion of the vulnerable capacity is 


* 
Ae exp(-b, sk; .) =U /d, < Val Bae { 


and the feasible minimum length for the route is 


* * 
= Gli eae Pr 1/b. 
ind is 
* 
The value of u- can now be computed in closed form and the value of. the 
feasible minimum length for the route obtained. The equation for com- 
* 
puting u is 


Wy = V) eCXp RISEN vee fie (1/b; 3) 


* 
Se leek is less thanmonmsecquad te h(v_). then u- is greater than or 
equal to Vy: A lower bound on the feasible route length is thus. pro- 
vided: 


—— ae) 


> a eee ee ee 1/b ay 
153 ie we mes) ini Sb aa 


This bound can prove useful in the following manner: if any route has 
a minimum feasible length or lower bound which is greater than the least 
upper bound for the network, this route can be eliminated from further 


consideration. 
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6) If any feasible route length is less than the least upper bound 
for the network, this becomes the new network least upper bound. 

7) Only those routes with lower bounds less than or equal to the 
Jeast upper bound of the network need be considered further. 

8) To continue, it is now necessary to perform the summations over 
the arcs which have Ws, greater than u for those routes for which 
u" and the minimum feasible length have not yet been determined. Begin- 
ning with the route having the next larger lower bound, order the Me 


in increasing order, where 


=Wb_= the q-th largest ~~ an : 
Vo a? q gest value of Wag? a3 in that route; 
fa Wb, = max erie : 
T5J 
-_ WD, = © Ws 3D, . 


Summations are henceforth made over the single subscript. Successively 


‘compute h(va) for q = n,n-1,...,]1 for the route being considered until 


h(v,) >K > h(v _)) For seme ty” “1 <q <M . 


q 


9) At each iterative computation, if it is found that h(v,) K, 
a new improved lower bound is provided as this means another arc has 
been found for which no sorties will be allocated and that u > ihe 
This means the feasible route length will not be as great as if u = Va? 
that is, 
q-1 qe] * 


n 
iS - >See > ee se. So (7) 
ee V4] 


If the lower bound for a route exceeds the current least upper 


bound for the network at any time, this route may be dropoed from further 


consideration. 


ZO 


10) Once the value of q is found so that both htv,) > K and 


K > n(Va_y)s the improved lower bound is found as in (7), and the least 


upper bound for the route 1s 


n q- | q-| 
ee a ee OT 
i=q i=] Se ate a 


* 
(1/b,) > M . (8) 
The above relation is a result of the fact that K greater than h(ve_y) 
* 
implies that u is less than Yi 


WG The minimum feasible route length can now be derived for the 


route being considered: 


" n of a 
M = a U. hi y L. + U > (1/b.) 
1=q i=] i=] 
where 
* q-| 
1 = ie. exp{ -[K - hve I E (1/b. ) |e 


i=] 

12) The process is repeated from step 6 until all routes have been 
considered. The route with the minimum feasible route length over al] 
routes in the network is the one to allocate airstrikes to. Ties indi- 
cate that there are two or more minimum feasible cut sets in the primal. 

13) Having selected the route which produces the nie feasible 
distance through the network and the u" for that route, the remainder of 
the solution is available. Those arcs with Hie Sloe greater than u" will 
be attacked and the level of attack on arc (i,j) will be vee from 


equation (6). 


D. EXAMPLE 
The following is a hypothetical example of a transportation network. 
The first diagram in figure | is the primal network. Flow is from node 


1 to node 6. 
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Figure 1. An Example Network 
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There are three numbers associated with each arc: U and Dig: 


iJ 
These represent the upper and lower bounds and the vulnerability para- 


ig? 


meter, respectively, of each arc. The number of sorties available for 
the interdiction of the network is 40. 

In figure 2 the topological dual is formed by first constructing the 
artificial arc between the source and the sink. Asterisks are used to 


denote the nodes of the dual. 


2] 





Figure 2. Construction of the Dual 


Figure 3 shows the topological dual with the appropriate figures 
transferred from the primal. For convenience of reference the various 
routes between the source and sink are designated by the nodes over 


which they pass. 


Ry: 1,2,5,8,11,14 Re: 1,256,9,7,10,13,14 
Ry: 1,2,6,9513,14 Ry: 1,3,6,2,5,8,11,14 
Ra: 1,3,6,9,13,14 Ro: 1,3,6,9,7510,13,14 
Ry: 1,4,7,9,13,14 Ry: Vs4at so sOs soso nie 4 
Re: 14457510,12,14 
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Figure 3. The Topological Dual 


The minimum length normal route is provided by R, with a length of 


] 
2130. This becomes the first least upper bound for the network. 

In determining the r shortest routes using the lower bounds, the 
following order is found: ; 

R, pede 1620, Re: 1770, R,: 1790, Ro: 1920, Ra: 2250 
The lower bound on the next shortest route is from Ro; re is 2250 which 
exceeds the present least upper bound. For this reason it is not neces- 
sary to consider Re RO Re or Rg: 

Computations are begun on the routes in increasing order of lower 
bounds. When a feasible minimum length is found which is lower than the 


lower bound of the remaining routes, the optimum solution has been found. 


No other route can be reduced to this amount. 


The computations for R are done first. The ordered Ve for R, are: 


a 


0.1(100) = 10; v, = 0.07(130) =9.1:V,= 0.08(100) = 8; 


V woe 


1 2 


0.06(50) = Saave= ORO yo = 17. 


“4 5 
The function n(v_) is computed, and it found that the number of 

sorties available is less than h(v_) = 70.51. Similarly it is found 

that K = 40 is less than h(vy) = 40.15, but K is greater than h( V2) 

= 4.07. This indicates that arcs (1,2), (11,14) and (8,11) would be 

attacked while arcs (2,5) and (5,8) would have no sorties allocated. 

A first improved lower bound ‘is obtained by adding the sum of the upper 

bounds of those arcs which would not be attacked and the sum of the 


lower bounds of those arcs which would be attacked if this route were 


chosen. 


Be = 


bt. =O 


] 


5 3 3 
fo. ou + yc. ee el) De) 
i= isl | 7 
= 1800 + 3(36.79) = 1910.37 
From (8) an upper bound for Ry is 
5 3 S 
Uae UE ae cL. gt v4 25 (/be 
i=q i=l | | re 
= 1800 + 8(36.79).= 2094.32 


This becomes the new least upper bound for the network and replaces the 
figure 2130 previously derived. 
The minimum feasible length for R, is 


* 
[is + U 
1 " 
| | 


(1/b; ) 


lm Gs 
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* 
where u_ is the solution to h(u) = K, 


* 


uy = vq exp [ - (40 - 4.0719)/36.7857] = 3.0125 . 
Therefore 
M = 850+ 950 + 3.01(36.79) = 1910.74 . 


This ‘is less than the current least upper bound and therefore becomes 

the new least upper bound for the network. R,, may now be eliminated from 
further consideration as it has a current lower bound of 1920. Only R25 
Ry and Re remain for consideration. 


The computations for Ra result in 


LS. = Wea , 
‘U.B. = 1990.70 , 
* 
mn i= “ae. 
* 
Mw = o0eRn4 . 


Ry is therefore better than R.. Only R. and Re remain for consideration 
as possible best routes. Computations for Re result in finding an 
amproved lower bound of 1936 and Re can be eliminated. Similarly the 
amproved Jower bound for R. is 2159 and R. need not be considered further. 
Ry as therefore the best route. The non-integer solution for the number 


of sorties to allocate to each arc is 


yp = 12.13, 
Tb) = 
kg.g = 11.55 
Lan ON 
ky qaz 10-93 - 


The optimal integer solution is 


he = &. 
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Oe 


WW 
— 
NO 


Keg 


Kea] OO 


Kia = IE 


which results in a feasible length of 1902.21 . This is better than the 
non-integer solution of the closest contender, Ry , 


Returning to the primal, the cut set to be attacked is as shown in 


figure 4. 





Figure 4. Minimum Cut Set After Attack 
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V. DISCUSSION 


A. UHITENESS 

The procedure described ends in a finite number of steps as the 
routes are systematically eliminated by means of upper and lower bounds, 
after which the worst possibility is complete enumeration for the re- 
maining routes. 

A non-integer solution is produced. If the parameters and the 
resource constraint are such as to be on the steep slope portion of the 
curve for the objective function, the integer solution may be vastly 
different than the non-integer solution. As the improved bounds in the 
procedure are for the feasible non-integer solution, it may be that the 
integer solution will actually be outside the non-integer bounds. If, 
however, the problem is such as to be ona less steep portion of the 
curve, the integer and non-integer solutions may vary only slightly. 

In this case a simple round-off procedure might be used to get an integer 
solution without much difference with the correct optimal integer solution. 
It appears that the only method to derive the optimal integer solution is 


by investigation around the non-integer solution. 


B. ASSUMPTIONS 

The assumptions that the values of maximum capacity, vulnerable 
proportions, lower bounds and vulnerability factors are known and deter- 
ministic is difficult to justify. Experimental testing such as that 
suggested by Holliday [13] might be used to determine whether the exponen- 
tial distribution could validly be used as a damage function. It would 


then be necessary to determine the probability of a target hit and the 


Sie 


proportion of capacity reduction done by a sequence of target hits under 
various conditions of road types, concealment, weather and other charac- 


teristics as might be judged important. 


C. INTERPRETATION OF THE LAGRANGE MULTIPLIER 

The Lagrange multiplier u has an interesting interpretation and 
consequent use. It is the differential rate of damage for an additional 
sortie. In this sense, it can be interpreted as the implicit marginal 
price which a commander 1s placing on the expected destruction of one 
more unit of network capacity when he limits the number of sorties avail- 
able for network interdiction. Another manner in which it might be used 
is when a commander decides he will] allocate strikes up to the point 
where he starts getting less than u units of capacity reduction for each 
Strike. In this sense he is making a trade-off with the other alterna- 
tives available to him. He allocates aircraft to tactical interdiction 
up to the point where the marginal return is in some sense the same as 
that of an additional aircraft allocated to any other mission. 

The solution to such a problem can be derived by using the given 
value of u in the formulae as the criterion for attack. Those arcs with 
Waa greater than u would be attacked at the level - where kegal'S 
determined from equation (6) using the given value of u as u 
D. ANOTHER JUSTIFICATION FOR USE OF THE EXPONENTIAL FUNCTION 

Koopman [9] gives an interpretation of optimum destructiveness and 
makes a comparison with search theory wherein he gives some justification 
for the use of the exponential as a damage function. He points out the 
problem of optimum distribution of gunfire, bombs or other destructive 
missiles projected into an area. His presentation will be paraphrased 
in terms of optimum allocation between arcs of a route. 
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Two arcs with maximum lengths (capacities in the primal) of U, and 
Un have proportional vulnerable lengths of Wo /U, and WA/U, respectively. 
Because of the type of bombing used, defenses, overhead cover or what- 
ever, the bombing will be sufficiently accurate to hit the particular 
arc but not to hit the vulnerable capacity in the area except by chance. 
The vulnerability of all capacity is assumed to be the same on a partic- 
ular arc. The result of n random hits on the target area W.,1=1,2, is 
assumed to reduce the vulnerable capacity by a proportion 1 - . where 
Ce Ts less than 1. Thus after n hits on the vulnerable capacity the 
total capacity of the arc would be reduced to U. - W. + ce W, , the 

n 


reduction being W. =a, W. and the proportionate reduction 


(W, 


- cy Ws), 

If there are N bombs available and no cost difference exists between 
dropping bombs on one arc than the other, then the decision as to how to 
al.locate bombs between the two arcs is properly made on the basis of 
maximizing the expected damage to be inflicted. 


If N, bombs are used on arc | and Ny bombs on arc 2, 


] 


N, +N, = N~ and N 


1 2 


The probability that one bomb dropped on the first arc will hit the 
target is W,/U; ,» assuming that the vulnerable capacity is uniformly 
distributed. (It is also possible to use as the probability of hit, p> 
instead of the ratio W./U..) If N, bombs are dropped on arc 1, the 


probability that R, of them will] hit the target is (by the binomial 


1 
distribution) equal to 


IM\/y 


\ R | | 
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The proportional damage done by the R hits being | - cu , the expected 


proportional damage produced by Ny bombs dropped on arc 1 is given by 


N 


\ / R \ 

i Ny Ww \" | ‘oy Ry ) 

/ ae mwa = ey = - exp{-b,N. } 
R,=0 oF / 


meres) = = as ill: Wi - C1 )/U, J : 


Thus the expected proportional damage is 1 - exp(-b,N.) . The 
expected reduction in capacity over the route when N, bombs are dropped 


on arc | and No bombs are dropped on arc 2 is 
WO - exp (-b)N, )] + WoL - exp (-bNo) J (9) 
The resulting expected capacity is 


U, - W [1 - exp (-b,N, ) + U, - W, exp (-b,N,) J (10) 


or equivalently 


L; + Wy exp (-b,N, ) + Ly + Ws exp (-b,N.) (11) 


Similar results are obtained by either maximizing (9) or minimizing 
either (10) or (11) subject to the restrictions that Ny + N, = N and 
N 


ae NN, > Os. 


1 2 


The bomb allocation problem is therefore to find the bombing strategy 
which gives the minimum expected value of the minimum feasible cut set 


through the flow network. 


E. APPLICATIONS 
Due to the particular characteristics of the transportation networks 


used in most of the land masses in Southeast Asia, the proposed procedure 
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might be quite appropriate for use in the interdiction problem in such 
areas. The fact that large portions of roadways in that area of the 
world are under a jungle canopy or can be easily provided with artificial 
concealment makes it difficult to obtain accurate bomb hits. This situ- 


ation lends itself to a probabilistic treatment. 


F. RECOMMENDATIONS FOR FURTHER STUDY 

An analysis of actual strike data is appropriate before the exponen- 
tial damage function can be considered for real-world use. Research such 
as that done by Holliday for ARPA [2,13] might provide such data. It 
would be desirable to determine hit probabilities and the expected 
reduction in capacity resulting from a sequence of hits by various weapon 
systems. 

The exponential function is particularly easy to use because it 
provides numerous simplifications to the computational problem. Gibbs! 
lemma however applies to any continuous differential function with 
decreasing marginal returns. The data might suggest that the exponential 
is inappropriate but that some other function having such marginal 
returns can be used. 

On the surface there appears to be some optimal criterion for the 
rounding-off of the non-integer solution. The branch-and-bound method 
does not seem to lend itself to a simple application in this case. 

As with many other treatments involving allocation, this paper was 
unable to deal with the cost of lost aircraft. In a way, this was taken 
care of in the vulnerability parameters of the various arcs. If the 
probability of getting through the defenses is 3/4 and the probability of 


a target hit is 1/3 for those aircraft getting through the defenses, then 


3/ 


the probability of a target hit by a sortie allocated can be expressed 


in terms of conditional probability as 


Wok 3 a = 1/4 
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VI. SUMMARY 


A procedure has been presented for computing the optimum allocation 
of aircraft sorties toward the tactical interdiction of a transportation 
network when the flow of supplies is restricted by the capacity of the 
network. A specific damage function is assumed which requires knowledge 
or estimates of the upper and lower bounds of capacities as well as the 
vulnerability of each arc to the weapon system involved. 

The method utilizes the notion of the topological dual and finds the 
minimum feasible route through the dual. The total length of this route 
corresponds to the minimum value to which the capacity of the network 
can be reduced given the number of sorties available. 

The alternative routes through the dual are first screened systemati- 
cally by means of bounding the minimum feasible length. The stepwise 
procedure eliminates those routes from consideration which cannot meet 
the first least upper bound. The Lagrange multiplier is used as a thresh- 
hold criterion based on the product of the parameters nee and Ds to 
determine whether an arc should be hit. The number of sorties to allocate 
to the arc is then determined. The routes are considered in increasing 
order of sums of lower bounds. Once a feasible length or least upper 
bound is found less than the next lower bound, all other routes may be 
eliminated from consideration. The route with the minimum feasible length 
is the best route and corresponds to the minimum cut set in the primal. 

In order to begin the stepwise procedure the inputs required are the 
upper and lower bounds on arc capacity, the vulnerability of each arc, and 
the total number of sorties available for interdiction of the network. The 
procedure determines which arcs to hit and at what level to achieve the 


minimum feasible capacity. 
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